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The dineutron correlations and the crossover from superfluidity of neutron Cooper pairs in the 
'So pairing channel to Bose-Einstein condensation (BEC) of dineutron pairs in both symmetric 
and neutron matter are studied within the relativistic Hartree-Bogoliubov theory, with the effec- 
£/■} [ tive interaction PK1 of the relativistic mean-field approach in the particle-hole channel and the 

finite-range Gogny force in the particle-particle channel. The influence of the pairing strength on 
the behaviors of dineutron correlations is investigated. It is found that the neutron pairing gaps 
at the Fermi surface from three adopted Gogny interactions are smaller at low densities than 
the one from the bare nucleon-nucleon interaction Bonn-B potential. From the normal (anoma- 
lous) density distribution functions and the density correlation function, it is confirmed that a 
true dineutron BEC state does not appear in nuclear matter. In the cases of the Gogny interac- 
tions, the most BEC-like state may appear when the neutron Fermi momentum kf„ ~ 0.3 frrT 1 . 
Moreover, based on the newly developed criterion for several characteristic quantities within the 
relativistic framework, the BCS-BEC crossover is supposed to realize in a revised density region 
£>. | with kf n € [0.15, 0.63] frrT 1 in nuclear matter. 
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_f~. , Pairing correlations are significant phenomena in the fermion systems. For the neutron- 

neutron pairing, the correlation is supposed to be important in deciding the properties of low- 
density nuclear matter. One of the evidence comes from the fact that the bare neutron-neutron 
interaction in the 1 Sq channel leads to a virtual state around zero energy characterized by a large 
negative scattering length a « -18.5 + 0.4 fm [1], implying a very strong attraction between 
two neutrons in the spin singlet state. Furthermore, theoretical predictions suggest that around 
1/10 of the nuclear saturation density po, the So pairing gap may take a remarkably larger 
value than that around po [2]. In addition, the strong dineutron correlations are also supported 
by the enhancement of two-neutron transfer cross sections in several heavy nuclei ||3j]. In the 
weakly bound neutron-rich nuclei, the couplings to the continuum could strengthen the dineutron 
correlations, which play a crucial role in unstable nucleus and the formation of the halo D4J,|5|]. 
Recently, a ground state dineutron decay was observed for the first time in 16 Be with a small 
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emission angle between the two neutrons, indicating the structure of dineutron clusters inside 
neutron-rich nuclei [6]. 

The BCS-BEC crossover is one of the hottest issues recently in the study of the pairing cor- 
relations in the fermion systems, which is described as a smooth and adiabatic evolution of the 
pairing from the weakly coupled Bardeen-Cooper-Schrieffer (BCS) type to the strongly corre- 
lated Bose-Einstein condensation (BEC) state with increasing pairing strength (17l|8[]. For the 
neutron-proton pairing, the BCS-BEC crossover has been investigated in the 3 S \- 3 D\ channel, 
in which the st rong spa tial correlation and the BEC of the deuterons may occur at low densities 
H OH [U EH 013 EI EH Eztl- On the other side, the progress in studying the dineutron 
correlations in weakly bound nuclei has stimulated a lot of interests in searching for possible 
BCS-BEC crossover of neutron-neutron pairing flH HH HI |H HHH Q . It has been 
revealed that the dineutron correlations get stronger as density decreases, and the BCS-BEC 
crossover of the neutron pairing may occur at low densities. The spatial structure of neutron 
Cooper pair wave function could evolve from BCS-type to BEC-type with decreasing density. In 
finite nuclei, the coexistence of BCS- and BEC-like spatial structures of neutron pairs has been 
exhibited in the halo nucleus n Li as well [26]. From the two-particle wave function, a strong 
dineutron correlation between the valence neutrons appears on the surface of the nucleus. 

As the relativistic mean-field (RMF) theory and its extension to the relativistic Hartree- 
Bogoliubov (RHB) theory had achieved lots of successes in the descriptions of both nuclear mat- 
ter and finite nuclei near or far from the stability line [5, 27], the dineutron BCS-BEC crossover 
in the 1 Sq pairing channel was studied within the RHB theory in nuclear matter, taking the bare 
nucleon-nucleon interaction Bonn-B in the particle-particle (pp) channel B24U25I1 . It is argued 
that there is no evidence for a true dineutron BEC state at any density, but some features of the 
BCS-BEC crossover are seen in the density regions, 0.05 fin < kf„ < 0.7 (0.75) fm" 1 , for the 
symmetric nuclear (pure neutron) matter. In addition, the reference values of several quantities 
characterizing the BCS-BEC crossover are obtained in the relativistic framework [25]. 

Since details of the nucleon-nucleon interaction are still not completely uncovered because 
of the insensitivity of the experimental data, for convenience, several effective nuclear interac- 
tions, such as the zero-range contact force or finite-range Gogny force [28], are widely used in 
the pp channel in most of investigations on nuclear pairing correlations. Besides, the effective 
interactions in the RMF theory are also used in the pp channel [29, 30], but one has to introduce 
an effective factor to this kind of pairing interactions in order to obtain reasonable values for 
the gap parameter 113011 - Because the pairing properties strongly depend on the pairing strength, 
it is expected that the BCS-BEC crossover phenomenon will be influenced by different selec- 
tions of the pairing interaction. It has been shown that the Gogny interactions predict pairing 
prop erties that are surprisingly close to those derived from the bare realistic nuclear interactions 
|3jJ,[32[]. Recently, a Gogny-Hartree-Fock-Bogoliubov nuclear mass model is presented, which 
makes great success in reproducing the nuclear masses with an accuracy comparable with the 
best mass formulas 03311 . Thus, it is of great interest to study the behavior of dineutron correla- 
tions and the BCS-BEC crossover in nuclear matter with several Gogny interactions used in the 
pp channel, and investigate the effects of different pairing forces on the results. 

It should be noticed that the bare nucleon-nucleon interaction should be corrected by the 
medium polarization effects (MPEs), including mainly self-energy contributions and vertex cor- 
rections, in extracting the effective interactions in the pp channel at various densities 134 l35l 
3a I37L l38l 13911 . Most of these studies predict a reduction of the pairing gap in neutron matter 



due to the MPEs. However, it has been suggested that the MPEs enhance the pairing interac- 
tion in finite nuclei because of the surface vibrations [40]. Encouraged by this apparent con- 
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tradiction between neutron matter and finite nuclei, a microscopic Brueckner calculation was 
performed on the MPEs in both symmetric and neutron matters and dramatic results were re- 
vealed [39]. The MPEs clearly increase the pairing correlations in symmetric nuclear matter 
(especially for neutron Fermi momentum kf„ < 0.7 far') but decrease them in neutron matter. 
Further studies claimed that a dineutron BEC state would be formed in symmetric nuclear mat- 
ter at kp n ~ 0.2 fm~' if the pairing interaction is screened II191 |20J. Moreover, the effects of 
medium polarization and possible ambiguity of pairing force on the dineutron correlations and 
the BCS-BEC crossover were investigated in the nuclear matter phenomenologically by taking 



the bare nucleon-nucleon interaction in the pp channel multiplied by an effective factor 125 1. 
It is shown that if the effective factor is larger than 1.10, a dineutron BEC state appears in the 
low-density limit, and if it is smaller than 0.85, the neutron Cooper pairs are found totally in the 
weak coupling BCS region. 

In this work, the dineutron correlations in the 'So channel for the nuclear matter will be stud- 
ied within the RHB theory, with three selected finite-range Gogny interactions, i.e., Dl [28], D1S 



114 111 and DIN 14211 . in the pp channel. The evolution of several characteristic quantities for the 
dineutron correlations, including the normal and anomalous density distribution functions as well 
as the density correlation function of the neutron Cooper pairs, will be discussed in comparison 
with the results with Bonn-B potential [24], then the BCS-BEC crossover will be investigated 
with the criteria based on the relativistic calculations N25I1 . As have been elucidated in the previ- 
ous study [25], the MPEs that go beyond the mean-field approximation are not considered in this 
work. Attention will be paid to the comparison among various versions of the pairing interaction. 
In the RHB theory, meson fields are considered dynamically beyond the mean-field theory to 
provide the pairing field through the anomalous Green's functions [29]. For the infinite nuclear 
matter, the Dirac-Hartree-Fock-Bogoliubov equation reduces to the usual BCS equation. For the 
'So channel, the pairing gap function A(p) is 

m ^^[ V ^ P) W dK (1) 

where v pp (k,p) is the matrix element of nucleon-nucleon interaction in the momentum space 
for the 'So pairing channel, and E^ is the quasi-particle energy. The corresponding normal and 
anomalous density distribution function are defined as 
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Kk = w (2) 



with the single-particle energy e* and the chemical potential p obtained from the standard RMF 
approach. For nuclear matter with given baryonic density p\, and isospin asymmetry f = (p„ - 
Pp)lpbi me g a P function can be solved by a self-consistent iteration method with no-sea approx- 
imation. For more theoretical details we refer the reader to Refs 151 124 12511 . 



The finite-range Gogny interactions are utilized in the pp channel, which are obtained by 
combining two Gaussian functions describing a long-range attraction and a short-range repulsion 
respectively. For the 'So pairing channel, the matrix element v pp (k,p) is just the average of 
v P p(k,p) over the angle between the vectors k and p, which has the form 

4tt 3 / 2 ^ _ / it .., ,.\ . . .1 
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Figure 1 : Neutron pairing gap at the Fermi surface A(kf „ ) as a function of the neutron Fermi momentum kf „ in symmetric 
nuclear matter (f = 0, left panel) and pure neutron matter (f = 1, right panel). The effective interaction PK1 |43] is used 
for the mean-field calculation in the ph channel, and the Gogny interactions Dl (dashed line), D1S (dashed-dotted line) 
and DIN (dotted line) are used in the pp channel in comparison with Bonn-B potential I24I (grey solid line). 



with the parameters 14 and C m = p,„(W m - B m - H m + M m ) taken from the corresponding Gogny 
interaction IJ28J l4ll 14211 . It is noteworthy that only the density-independent part of the Gogny 
interaction contributes to the pairing in the So channel. 

In the following studies, the Gogny interactions Dl, D1S and DIN will be adopted for 
v pp (k,p). For the mean-field calculation in the particle-hole {ph) channel, the effective inter- 



action PK1 1 43] is used, since the results do not depend sensitively on various other parameter 
sets 1 24]. For comparison, the results given by the pairing force with bare nucleon-nucleon 
interaction Bonn-B [24] will be discussed again as well. 

To investigate the properties of the neutron Cooper pairs, it is interesting to look into the 
neutron pairing gap A(p) first. One of the most important properties of the neutron pairing 
gap is its value at the Fermi surface Af„ = A(kf n ). In Fig. [1] Af„ is shown as a function of 
the Fermi momentum kf n for different pairing interactions in both symmetric nuclear matter 
and pure neutron matter. It is found that the neutron pairing gap Af„ is strongly dependent on 
the Fermi momentum, or equivalently, the nuclear matter density. Af n increases as the Fermi 
momentum (or density) goes down, reaches a maximum at kf n « 0.8 frrT 1 in symmetric nuclear 
matter or kf„ « 0.9 firT 1 in pure neutron matter, and then rapidly drops to zero. A systematical 
enhancement of about 0.3 MeV for Af„ around kf„ - 0.8 fm -1 is revealed in pure neutron matter 
compared with those in symmetric nuclear matter for all of the adopted pairing interactions. 

In comparison with the Bonn-B potential, it is seen that the Gogny interactions especially Dl 
have larger pairing gap when approaching to the saturation density. At kf„ = 0.8 fm" , among 
the adopted pairing interactions, Dl gives the maximum values of Af„, namely, 3.13 MeV in 
symmetric nuclear matter and 3.40 MeV in pure neutron matter. As the density decreases, the 
values of Af„ from Dl, D1S and DIN become first consistent with and then smaller than the 
one from Bonn-B potential gradually. As illustrated in a recent work [25], the pairing strength 
required for appearance of a dineutron BEC state in the low-density limit must be stronger than 
1 . 1 times of Bonn-B potential, and the corresponding pairing gap Af„ is 4. 12 MeV around kf n - 
0.8 fm -1 . Therefore, it is expected here again with the Gogny pairing interactions that a true 
dineutron BEC state cannot occur at any density in nuclear matter. 

The normal and anomalous density distribution functions, pk and Kk, provide us valuable 
informations on the mean field and the pairing field. In Fig. [2] neutron normal and anomalous 




Figure 2: Neutron normal and anomalous density distribution function, pt and Kk, as a function of the ratio of the neutron 
momentum to the Fermi momentum k/kf „ at several neutron Fermi momenta kp n in pure neutron matter. The effective 
interaction PK1 ["43 ] is used for the mean-field calculation in the ph channel and the Gogny interaction Dl is used in the 
pp channel. 



density distribution functions are plotted as a function of k/kf„ at several neutron Fermi momenta 
kf„ in pure neutron matter with the Gogny pairing interaction Dl. When kf„ = 1.4 fm , the 
normal density distribution function pk is represented nearly as a step function at k/kf„ = 1.0, and 
the neutron Cooper pair could be treated as a BCS-like pair. With decreasing Fermi momentum, 
the momentum distribution of pk evolves smoothly and deviates from the step function gradually. 
After reaching a maximum deviation from the step function at kf„ - 0.2 far 1 and 0.4 far 1 , the 
momentum distribution of pk approaches to the step function again at dilute density. 

The anomalous density distribution function Kk, known as the Cooper pair wave function 
in momentum space, presents a similar density-dependent behavior to pk as well. When kp n = 
1 .4 firT 1 , it is seen that the pairing exists almost on the Fermi surface with a well-developed peak 
of K/t. As the density decreases, the momentum distribution of Kk is expanded to both lower and 
higher regions than the Fermi momentum. After generating a maximum momentum dispersion 
of ^ at kf„ = 0.2 firT 1 and 0.4 far 1 , the momentum distribution of k> is narrowed again in 
the vicinity of the Fermi momentum at dilute area. Therefore, from the evolution of pk and Kk 
with the Fermi momentum, the most dineutron BEC-like state may appear at the density region of 
kp n G [0.2, 0.4] far 1 in the case of the Gogny pairing interaction Dl. Similar results are obtained 
with the other two Gogny pairing interactions D1S and DIN, and for symmetric nuclear matter 
as well. 

As a quantitative study, the density correlation function developed recently for ultracold 
atomic gases has been suggested as a clear distinction between the BCS and BEC limits, which 
gives prominence to the relative strength between the mean field and the pairing field [20, 24, 44] . 
It is expressed as 

D(q) = I K (q) - I p (q). (4) 
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At zero-momentum transfer q — 0, the normal and anomalous density contributions I p and l K 
respectively read 

1 C°° 

7,(0) = -=- p^dfc, (5) 

"P» Jo 

1 f°° 
4(0) = — 4k 2 dk, (6) 

7rp» Jo 

and they satisfy the sum rule / p (0) + 4(0) = 1 . The sign change of the density correlation function 



at zero-momentum transfer D(0) could be taken as a criterion of the BCS-BEC crossover [44], 
i.e., D(Q) < means a BCS-type pairing and D(Q) > represents a dineutron BEC state. 

In Figj3] the density correlation functions at zero-momentum transfer D(Q) and its normal 
and anomalous density contributions, 7 p (0) and I K {0), are shown as a function of kf„ in nuclear 
matter. In three cases of Gogny interactions, as the density decreases, it is seen that 7 p (0) first 
drops down from 1.0, reaches a minimum at kp„ ~0.3 firT 1 , and then goes up again, while I K (Q) 
gives opposite trend and has no contribution when approaching to the saturation density. Thus, a 
peak for the density correlation function D(0) is found around kf„ = 0.3 fm _1 for the Gogny in- 
teractions, larger than the value of kf„ - 0.2 fin" 1 in the case of the Bonn-B potential. Therefore, 
the most BEC-like state may appear around kf„ - 0.3 firT 1 in the calculations with the Gogny 
interactions, which is in agreement with the former prediction by the normal and anomalous den- 
sity distribution functions, but distinguished from the Bonn-B case [24]. Because the anomalous 
density contribution I K (Q) is always smaller than the normal one 7 p (0), the density correlation 
function D(0) never changes sign at all for all of adopted pairing interactions. According to the 
criterion mentioned above, dineutrons are not in BEC state but just in the transition region be- 
tween BCS and BEC regimes. Compared with the result of Bonn-B potential, the values of D(0) 
are smaller in the cases of three Gogny interactions so that the corresponding dineutron behavior 
is farther away from the BEC limit. A further verification is performed here by the quasi-particle 
excitation spectrum with the similar method discussed in Ref. (24J, and the consistent conclu- 
sion that there is no evidence for the appearance of a true dineutron BEC state at any density is 
gained. 

As discussed above, the Gogny pairing interactions have lead to several recognizable differ- 
ences from the Bonn-B potential for the neutron pairing properties in nuclear matter. Hence, 
it is desirable to investigate the influence of the pairing interaction on features of the BCS- 
BEC crossover. Lots of characteristic quantities have been introduced to describe the BCS-BEC 
crossover quantitatively 1U8J ll9J |24J, |45|[. The probability P{r) for the pair partners within a rel- 
ative distance r is such a quantity to clarify the spatial correlation of the neutron Cooper pairs, 
which is defined as 

P(f) = f \V pair (r')\ 2 r' 2 dr', (7) 

Jo 

where v F pa ir('") is the wave function of the neutron Cooper pairs in coordinate space, 



TW(r) = -^-3 [ K k e ikr dk, (8) 

V-ny J 



with the constant C determined from the normalization condition. 

In Fig|U the probability P{d„), i.e., P(r) for the pair partners within the average inter-neutron 
distance d„ = p„ is shown as a function of kf„ with the different Gogny pairing interactions 
in both symmetric nuclear matter and pure neutron matter. It is revealed that the values of P(d n ) 
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Figure 3: Zero-momentum transfer density correlation function D(q = 0) and its normal and anomalous density contribu- 
tions, /p(0) and 7,(0), as a function of the neutron Fermi momentum kf„ in symmetric nuclear matter (f = 0, left panels) 
and pure neutron matter (£ = 1, right panels). The effective interaction PK1 |43] is used for the mean-field calculation in 
the ph channel, and the Gogny interactions Dl (dashed line), D1S (dashed-dotted line) and DIN (dotted line) are used in 
the pp channel in comparison with Bonn-B potential |24] (grey solid line). 
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Figure 4: Probability P(r) for the partner neutrons correlated within the average inter-neutron distance r = d„ as a 
function of the neutron Fermi momentum kf„ in symmetric nuclear matter (f = 0, left panel) and pure neutron matter 
(£ = 1, right panel). The effective interaction PK1 |43] is used for the mean-field calculation in the ph channel, and the 
Gogny interactions Dl (dashed line), D1S (dashed-dotted line) and DIN (dotted line) are used in the pp channel. The 
referred BCS-BEC crossover regions [25] (grey grid) are given as well. 



increase monotonically with decreasing density, and approach to 1.0 at dilute area. From the 
newly developed criterion based on the RHB theory, it is suggested that the BCS-BEC crossover 
is determined by P(d n ) > 0.81 in symmetric nuclear matter and P(d„) > 0.82 in pure neutron 
matter respectively 02511 . This limits the Fermi momentum for the realization of the BCS-BEC 
crossover in kf„ < 0.75 (0.78) frrT 1 (parameter set Dl), kf„ < 0.70 (0.75) firT 1 (parameter set 
D1S), and kf„ < 0.65 (0.70) firT 1 (parameter set DIN) in symmetric nuclear (pure neutron) 
matter, as shown in the columns marked with UL* in Table Q] 
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Figure 5: The ratio Af„/e/r„ between the neutron pairing gap at the Fermi surface Af„ and the neutron Fermi kinetic 
energy e? „ as a function of the neutron Fermi momentum kp„ in symmetric nuclear matter (f = 0, left panel) and pure 
neutron matter (f = 1, right panel). The effective interaction PK1 [43] is used for the mean-field calculation in the ph 
channel, and the Gogny interactions Dl (dashed line), D1S (dashed-dotted line) and DIN (dotted line) are used in the pp 
channel. The referred BCS-BEC crossover regions (grey grid) and the unitary limits (dotted line) 12511 are given as well. 



The BCS-BEC crossover could be explored as well by the ratio Ap n jep n between the neutron 
pairing gap at the Fermi surface Af„ and the neutron Fermi kinetic energy ef n = e(kf„). Here the 
neutron kinetic energy e(k) is given by 



e(k) = V/fc 2 + M* 2 - M*, (9) 

where the scalar mass M* can be obtained from the standard RMF approach [5]. If the ratio 
is large enough, the neutron is expected to be in the strong coupling regime. From the criterion 
based on the RHB theory, the BCS-BEC crossover region is estimated as 0.26 < Af„/ef„ < 1.50 
in symmetric nuclear matter and 0.28 < Af„/ef„ < 1 .52 in pure neutron matter, while the unitary 
limit is given by Af„/ef„ « 0.82 (0.83) in symmetric nuclear (pure neutron) matter. One should 
notice that the values of the criterion given here are signally different from those obtained in 
the regularized contact interaction model of non-relativistic framework 11181 Il9ll . Therefore, the 
density region for the realization of BCS-BEC crossover is supposed to change more or less 
compared to the previous results 112411 . even determined differently for symmetric nuclear matter 
and pure neutron matter. 

In Fig(5J the ratios Af„/ef„ are plotted as a function of the neutron Fermi momentum in 
nuclear matter based on three adopted Gogny interactions. With decreasing Fermi momentum, 
the values grow up first and go into BCS-BEC crossover region. It is found that the curves never 
intersect with the line of the unitary limit, i.e., the dineutron cannot perform a BEC state. At very 
low densities, the curves bend downward and return back to the BCS region again. Estimated 
from the criterion of Af„/ef„ referred above, the BCS-BEC crossover is marked in the density 
region with 0.08 (0.10) fin -1 < k Fn < 0.72 (0.73) fm" 1 (parameter set Dl), 0.10 (0.12) for 1 < 



Table 1: The lower (LL) and upper limits (UL) of the neutron Fermi momentum kf n for the realization of BCS-BEC 
crossover determined by the referred criteria of the ratio Af„/e f „, and the upper limits (UL*) determined by the referred 
criteria of the probability P(d n ) 1 25], with different pairing interactions in symmetric nuclear matter (f = 0) and pure 
neutron matter (f = 1). The values have units of fnT 1 . 
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k Fn < 0.70 (0.70) frrT 1 (parameter set D1S), and 0.14 (0.15) frrT 1 < k Fn < 0.63 (0.63) frrT 1 
(parameter set DIN) respectively for symmetric nuclear (pure neutron) matter. The results are 
summarized in Table [1] where the values determined from the newly developed criteria are also 
listed for the bare nucleon-nucleon interaction Bonn-B. 

From TableQ"! it is shown that the density region for the realization of the BCS-BEC crossover 
is affected clearly by the pairing force. In fact, the Gogny interaction DIN gives the weakest 
pairing gap at low densities as revealed in Fig. [T] so it generates a narrowest range of density 
for the BCS-BEC crossover, namely, 0.14(0.15) firT 1 < kf„ < 0.63 frrT 1 for symmetric nuclear 
(pure neutron) matter. Therefore, after combining the above analyses from the probability P(d n ) 
and the ratio Af„/ef„, the BCS-BEC crossover is supposed to appear most likely in the density 
region with kf„ e [0.15,0.63] firT 1 in nuclear matter, supported by not only three effective 
Gogny interactions but also bare nucleon-nucleon interaction Bonn-B. Additionally, it should 
be announced that if one utilizes the criterion obtained from the regularized contact interaction 
model in the non-relativistic framework 11181 Il9fl . the BCS-BEC crossover region will turn out 



to be kf„ € [0.11,0.72] firT 1 , a little broader than the one from newly developed criteria just 
mentioned. 

As a useful supplement, it is worthwhile to have a discussion about the role of the scattering 
phase shifts in the pairing properties. Several commonly used bare or effective nuclear interac- 
tions in the pp channel are deduced virtually from fitting the measured data of scattering phase 
shifts 1 191 I46I1 . Thus, corresponding pairing properties obtained are in agreement with those 
from quantum Monte Carlo calculations more or less, such as in Ref. 113811 . Several efforts have 
been made to give a direct analysis of the dineutron correlations in the low-density region from 
the scattering data [18, 19, 25], where the scattering length a, defined in terms of the T-matrix 
for the scattering in the free space, was used. According to the definition, the negative value 
of the scattering length a indicates an unbound state of neutron Cooper pairs while the positive 
value represents a bound state. It has been proposed to define the boundaries of the BCS-BEC 
crossover using a regularized gap equation approach fl7j [l8l |25j|45J], which is related with the 



scattering length a. It was shown that the properties of pairing correlations can be uniquely con- 
trolled by a dimensionless parameter l/(kf n a) which can give the evolution from BCS to BEC. 
The dimensionless parameter l/(kf n a) <K -1 corresponds to the weak coupling BCS regime, 
while l/(kf„a) » 1 is related to the strong correlated BEC regime. The boundaries characteriz- 
ing the BCS-BEC crossover can be approximately determined by 1 j(kf n d) — + 1 , and the unitary 
limit is defined as l/(kf n a) - 0, which is the midpoint of the BCS-BEC crossover. With this 
method, the BCS-BEC crossover region of the neutron pairing was determined within different 
many -body theories 1181 Il9l 12511 . Similar investigation is performed here with three selected 

9 



Gogny pairing interactions as well, and the results show good consistence with above analysis of 
the ratio A F „/e F „. 

However, one should be aware that some approximations have been introduced into the cal- 
culations of the regularized gap equation from the current RHB theory. The regularized gap 
equation approach is genetically applied to dilute systems, where the interaction matrix elements 
v pp (k, p) can be approximately treated as constant vo, while in the RHB calculations v pp (k, p) and 
correspondingly the pairing gap A(k) are found to be momentum dependent. Thus, such regular- 
ized approach gives poor knowledge on pairing properties at high densities. On the other hand, at 
very low densities the mean-field approach seems to be inefficient and some new features of nu- 



clear matter such as clustering should be considered [47 48]. More sophisticated approaches for 
the low-density nuclear matter is from the quantum Monte Carlo techniques 11221 l23l l38l l49l dOII 
or from the virial expansion method 115 ll 15211 . which are the best suited tools for treating strongly 
correlated systems thus could give reliable discussion on the large scattering length physics. By 
comparing the curves of the normal density distribution function p^ in Fig. [2] with those from 
quantum Monte Carlo calculations as shown in Ref. B49I1 . one find the results in Fig. [2] are lo- 
cated between the case of the step function and of the unitary limit, which confirm again that 
dineutrons at low densities are not in BEC state but just in the transition region between BCS 
and BEC regimes. In the virial expansion method, the second virial coefficient is related to the 
nucleon-nucleon scattering phase shifts so that the equation of state of low-density nuclear mat- 
ter can be derived directly from the measured scattering data. Therefore, it is deserved to take a 
look at the pairing properties with such kind of method. 

In conclusion, the dineutron correlations and the BCS-BEC crossover phenomenon for nu- 
clear matter in the ' S o channel have been investigated based on the relativistic Hartree-Bogoliubov 
theory with the effective RMF interaction PK1 in the ph channel and the finite-range Gogny force 
in the pp channel. The influence of the pairing strength on the behaviors of dineutron correla- 
tions is performed by comparing the results given by three selected Gogny interactions, i.e., Dl, 
D1S and DIN, and bare nucleon-nucleon interaction Bonn-B as well. It is found that the neutron 
pairing gaps at the Fermi surface Ap„ from three adopted Gogny interactions are smaller than 
those from the Bonn-B potential at low nuclear densities. From the normal and anomalous den- 
sity distribution functions, p^ and k^, and the density correlation function D(0), it is confirmed 
that a true dineutron BEC state does not occur at any density in nuclear matter, and neutron pair- 
ing is just in the region of BCS-BEC crossover at low densities, in agreement with our previous 



studies 124 12511 . However, from the calculations with the Gogny interactions, the most BEC-like 
state may appear at kf„ ~ 0.3 fm , which is larger than the value of kf„ ~ 0.2 fm" 1 in the case 
of the Bonn-B potential. Moreover, based on the newly developed criteria for the characteristic 



quantities B25I1 . namely, the probability P(d„) and the ratio Af„/ef„, the BCS-BEC crossover is 
marked in the density region with kf„ e [0.15, 0.63] firT 1 in nuclear matter, showing a narrower 
and more believable range than kf„ e [0.1 1, 0.72] firT 1 determined from the regularized contact 
interaction model in the non-relativistic framework and kf n e [0.05, 0.70] fm -1 in our previous 
study as well 10]. 
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